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Abstract 



D , The spin of supersymmetric particles can be determined at e + e colliders unambiguously. This is 

demonstrated for a characteristic set of non-colored supersymmetric particles - smuons, selectrons, 
and charginos/neutralinos. The analysis is based on the threshold behavior of the excitation curves 
' for pair production in e + e~ collisions, the angular distribution in the production process and decay 

angular distributions. In the first step we present the observables in the hclicity formalism for 
the supersymmetric particles. Subsequently we confront the results with corresponding analyses of 
Kaluza-Klein particles in theories of universal extra space dimensions which behave distinctly different 
from supersymmetric theories. It is shown in the third step that a set of observables can be designed 
which signal the spin of supersymmetric particles unambiguously without any model assumptions. 
Finally in the fourth step it is demonstrated that the determination of the spin of supersymmetric 
particles can be performed experimentally in practice at an e + e~ collider. 



1 INTRODUCTION 



The spin is one of the characteristics of all particles and it must be determined experimentally for any 
new species. Compelling arguments have been forwarded which suggest the supersymmetric extension 
of the Standard Model [1,2]. In supersymmetric theories (SUSY) spin-1 gauge and spin-0 Higgs bosons 
are paired with spin-1/2 fermions, gauginos and higgsinos, which mix, in the non-colored sector, to 
form charginos and neutralinos. Analogously spin-1/2 leptons and quarks are paired with spin-0 scalar 
sleptons and squarks. This opens a wide area of necessary efforts to determine the nature of the new 
particles experimentally. 

Measuring the masses of the particles is not sufficient to unravel the nature of the particles and of 
the underlying theory. This point has been widely discussed by comparing supersymmetric theories with 
theories of universal extra space dimensions (UED) [3,4] in which the counterparts of the supersymmetric 
partners are Kaluza-Klein (KK) excitations of the standard particles. When supersymmetric squarks 
are produced [5] at LHC, they may cascade down [6] to standard particles in the chain q — > qx® — * 
qli — > qiixl, which generates the observable final state qll. However, an analogous cascade can be 
realized in theories of universal extra space dimensions, starting from a KK excitation q\ of a quark, 
q\ — ► qZ\ — > qti\ — > qilji [7]. The origin of the observed chain particles, supersymmetry or extra space 
dimensions, can clearly be unraveled by measuring the spins of the intermediate cascade particles. 

Spin measurements of supersymmetric particles are difficult at LHC [7-9]. While the invariant 
mass distributions of the particles in decay cascades are characteristic for the spins of the intermediate 
particles involved, detector effects strongly reduce the signal in practice. 

In contrast, several techniques can be exploited to determine unambiguously the spin of particles 
produced pairwise in e + e~ collisions. These techniques have first been worked out theoretically for 
Higgs bosons, studied in the Higgs-strahlung process [10]; subsequent experimental simulations have 
proven these techniques to work in practice [11]. To conform with its scalar character, the polar angle 
distribution in smuon pair production has been investigated directly by reconstruction in Ref. [12] and 
reflected in their decay products in Ref. [13] at TeV and multi-TeV e + e~ colliders, respectively. 

A sequence of techniques, increasing in complexity, can be applied to determine the spin of particles 
in pair production 

e + e~ -> p, + pT, e + e~ and x + X~, X°X° (1-1) 

of sleptons, charginos and neutralinos in e + e~ collisions: 

(a) rise of the excitation curve near the threshold; 

(b) angular distribution in the production process; 

(c) angular distribution in decays of the polarized particles, 

eventually supplemented by 

(d) angular correlations between decay products of two particles. 
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While the second step (b) is already sufficient in the slepton sector, only the final state analysis is 
sufficient in general, including charginos/neutralinos, to determine the spin unambiguously. On the 
experimental side we follow the standard path. It will be shown in detailed simulations that the theo- 
retically predicted distributions in supersymmetric theories can be reconstructed after including initial 
and final state QED radiation, beamstrahlung and detector effects. Within the extended theoretical 
frame it is then proven that the assignment of the spin is unambiguous indeed. 

The report is organized as follows. In the subsequent Sections 2 to 4 we set up the theoretical basis 
for spin measurements of smuons, selectrons, and charginos/neutralinos. The technical frame we have 
chosen is the helicity formalism. We analyze which observables must be measured to determine the spin 
unambiguously. Moreover, simulations will assure us that the analyses of supersymmetric theories in 
e + e~ collisions can be performed experimentally. In the last Section 5 we briefly summarize the results. 
General formulae for the production cross sections of supersymmetric particles in collisions of polarized 
electrons and positrons are presented in an Appendix. 

2 SPIN OF SMUONS 

2.1 Smuon Production in e + e~ Collisions 

Smuons are the prototype for scalar particle pair production in e + e~ collisions [14, 15] mediated by 
the s-channel exchange of 7 and Z boson. Different lepton numbers prevent the flow of particles from 
the initial to the final state. For the sake of [experimental] simplicity we will restrict ourselves to the 
analysis of R-type smuons, 

e + e~ -> p&iiR (2.1) 

as these particles almost exclusively decay through the 2-particle channel p,^ —>■ [li^Xi with only one 
escaping invisible particle. The process is described by diagram (a) in Fig. 1. 




Figure 1: (a) s-channel 7 and Z exchange diagrams contributing to the production of all slepton pairs 
in e + e~ annihilation; and (b) t-channel neutralino exchange diagram contributing only to the production 
of selectron pairs in e + e~ collisions. 

The amplitude describing this production process can be expressed in terms of the generalized 
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electron charges 

Ql = 1 + (4/ - 1/2) c^D z ( S ) (2.2) 
Q R = l + t 2 w D z (s) (2.3) 

with s 2 ^ = sin 2 9w etc, 6w being the electroweak mixing angle, and the normalized Z propagator 
Dz{s) = s/[s — m?7 + imzTz], s denoting the squared center-of-mass energy [Dz is approximately 
real in the high energy limit s S> rn^]. The indices L and R in Eqs. (2.2) and (2.3) refer to left- and 
right-handedly polarized electrons [and oppositely polarized positrons], respectively. 

The total cross section and the distribution in the polar angle 9 between the flight direction and 
the e + e~ beam axis can be written in the form 1 

a[e+e-^Afc] = ^ [Ql + Ql] (2-4) 

--T^—R \ e+e ~ ->■ P-rPr] = Tsin 2 ^ (2.5) 

The coefficient /3 3 , with j3 = (1 — 4m 2 ifl /s) 1//2 denoting the velocity of the smuons, is the product of the 
phase space suppression factor 0, and the square of the P-wave suppression ~ near the threshold. 
The scalar smuon pair is produced in a P-wave to balance the spin 1 of the intermediate vector boson. 
Angular momentum conservation leads also to the sin 2 dependence of the differential cross section as 
forward production of spinless particles is forbidden. 

For asymptotic energies the cross section 

a — > -, for s — > 00 (2.6) 

follows the appropriate scaling law. 

The production of spin-0 particles in e + e~ annihilation is thus described by two characteristics: 

#1 threshold excitation ~ /3 3 (2.7) 
#2 angular distribution ~ sin 2 (2.8) 

The threshold excitation for smuons and the angular distribution are illustrated in Fig. 2(a) and (b), 
respectively. The (1r mass is chosen rrif LR = 300 GeV. In the following subsections it will be proven that 
the angular distribution is characteristic indeed for spinless particles and that it can be measured with 
great accuracy in e + e~ collisions. 



2.2 KK Excited States fif in UED 

The minimal UED version with one universal extra space dimension, which is compactified on the 
orbifold S1/Z2, generates a tower of spin-1/2 KK states over each L-type and R-type fermion in the 

1 The complete set of 1-loop radiative corrections, including the genuine SUSY corrections has been presented in Refs. [16, 
17]; see also Ref. [18]. 
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Figure 2: (a) The threshold excitation for smuons; and (b) the angular distribution in e + e~ — ► p^p^ 
[arbitrary normalization]. The predictions for smuons are compared with the corresponding observables 
for the first KK excitation in e + e~ — > MjiiMiji- 

Standard Model (SM) without generating additional zero modes [3]. Focussing on the R-type states 
in analogy to the previous subsection, the process for the production of a pair of the first muonic KK 
excitation, 



e e 



(2.c 



is described by the same diagram Fig. 1(a). fj,^ is a massive Dirac fermion carrying the electroweak 
charges of //^ but coupling only through vector currents to the electroweak gauge bosons 7 and Z. Thus 
the generalized charges are identical with Eqs. (2.3) and (2.2). It decays into the channel where 

71 is the lightest and stable excited gauge boson generally identified with the U(l) boson. 2 

The total cross section and the angular distribution are markedly different from the supersymmetric 
case: 
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(3-/? 2 ) 



[Ql + Q 



[l + cos 2 fl + (l-/3 2 )sin 2 9] 



(2.10) 
(2.11) 



8 (3 - /? 2 ) 

[Because of the vector character of the electroweak currents the distribution is symmetric in forward 
and backward direction.] The /iri pair is produced in an 5-wave with non-vanishing amplitude at the 
origin. The onset of the excitation curve is therefore suppressed only by the phase space factor ~ j3. The 
angular distribution is familiar from QED, being isotropic at threshold and evolving to the transverse- 



2 For the sake of simplicity we ignore the electroweak mixing of the KK excited gauge bosons and identify = Zi 
and Bi — 71 in the SU(2) and U(l) gauge boson sectors; the mixing is suppressed if the KK scale is much larger than the 
electroweak scale [4]. 



polarization term (1 + cos 2 9) asymptotically. 
For high energies the total cross section 



207ra 2 1 



for s -> oo (2.12) 



scales in the same way as the cross section but with a coefficient 4 times as large, as familiar from 

QED processes. 

Choosing for illustration a mass of 300 GeV, a value at the lower limit of the experimentally 
allowed range [3], the onset of the cross section and the angular distribution are displayed in Fig. 2(a) 
and (b). The center-of-mass energy is set to y/s = 1 TeV for ILC [14] and 3 TeV for CLIC [15] in 
the second figure. Both characteristics are markedly different from supersymmetric theories; see also 
Ref. [13]. In contrast to scalar smuons the onset of the excitation is vertical, proportional to the velocity 
(3 ~ [s — 4m 2 tJ?1 ] 1 / 2 as familiar for spin-1/2 particles. Also the angular distributions for scalar smuons 
and fermionic KK muon states are distinctly different. 

Already from the totally inclusive measurement of the production cross section no more confusion 
can arise between supersymmetric theories and theories of universal extra space dimensions. 



2.3 General Analysis 

Though the difference between the characteristics in the production of supersymmetric scalar particles 
and KK excited fermions can be exploited to rule out the false theory experimentally, we should explore 
nevertheless whether the conditions (2.7) and/or (2.8) are not only necessary but also sufficient to single 
out the scalar solution. 

The general analysis is most transparent if performed in the helicity formalism. In the process 

e + e - _ F J Xi f( 2 (2.13) 

for pointlike particles 3 F J and F with spin J and helicities Ai and A2 [either half-integer or integer], 
and mediated by s-channel 7 and Z exchange, the helicities of the electron and positron in the initial 
state are coupled to a vector with spin m = ±1 along the beam axis. The right-hand side of the diagram 
in Fig. 1(a) may then be interpreted as the decay of a virtual vector boson with polarization m = ±1 
to the F J F J pair. If the flight axis of the F J particles includes the angle 9 with the vector boson 
polarization axis [identical with the beam axis], the decay amplitude may be expressed in term of 
the helicity amplitude 7xiA 2 > 

{F( i Fl 2 -9\V m ) = VYsg v <L tW {9)T XlX2 with A' = Ai - A 2 (2.14) 

The angular dependence is in total encoded in the Wigner d function while the reduced helicity ampli- 
tudes T\ 1 \ 2 are independent of m. The value of A' is restricted to A' = and ±1. The totally inclusive 

3 Though only string theories are known to be consistent for interacting particles with J > 2, weakly interacting field 
theories can nevertheless be studied in approaches as formulated in Ref. [19]. 
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cross section can be expressed in the form 
ira 2 



(|?A,A-l| 2 + |T\,A+l| 2 ) + ^ l'^ XA l 



A A 

and the forward-backward symmetric part of the differential cross section analogously as 



da' 7T& 



dcos9 As 



P{Ql + Q\) 



1 + cos 2 9 



(I T a,a-i| 2 + |T A , A+1 | 2 ) + sin 2 6 Y, I^a| 2 



(2.15) 



(2.16) 



summed over the helicities of the outgoing F J and F J particles. To evaluate these expressions, the two 
cases in which F J is either fermionic or bosonic must be distinguished. 



(a) Fermionic spectrum F[J = 1/2, 3/2, . . .]: 

For pointlike theories in which the fermions carry electric and weak monopole charges the associated 
vector current includes the basic component [19] 



ft = 1p ai -a n l^ 



(2.17) 



built by the J = n + 1/2 spinor-tensor wave-function ip ai >— > Qn [20]. In analogy to the Dirac spin 1/2 
case, the current (2.17) can be decomposed, cf. Refs. [19,21], into an electric current and a magnetic 
current j™ as j M = j* + j™ with 4 



^ 2m' 

reducing to electric monopole and magnetic dipole currents in the non-relativistic limit. 



(2.18) 
(2.19) 



Both the electric current and the magnetic current j™ give rise to diagonal reduced helicity 
amplitudes through spin-O/D-wave and spin- 1/5- wave interactions, respectively; apart from overall 
coefficients, 



AA 



V2 

rrm _ T 



7 /3 2 \ (J + X) j-i/2 (J-A) j-i/2 

- Qx-l/2 (7) h~J~ Qx+1/2 (v 



2J 

( J + A ) o J ~ 1/2 h) - (J ~ A) o J - 1/2 h) 



(2.20) 
(2.21) 



Only the magnetic dipole current j™ generates non-diagonal reduced helicity amplitudes, 



A,A±1 



V(JTA)(J±A + 1) J-i/2 M 



(2.22) 



4 Demanding asymptotic unitarity at high energies, additional terms must be included in the basic Lagrangian [22] 
which alter the gyromagnetic ratio from g — J -1 universally to g = 2 [19], i.e. the coefficient in front of Eq. (2.19) may be 
adjusted accordingly. Universality of the value g = 2 beyond the Dirac theory is a well-known prediction of non-abelian 
gauge theories. 
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Here, 7 = y / s/2m is the Lorentz boost factor of the final-state particle and the energy-dependent 
functions (7) (|n| < N for integral N and n) are defined as 

N 2 N {N + n)\{N-n)\ _ _ A (27 2 ^o - 1) , . 

Qn(7) - ^ 2. • • W~+*»r 11 (1 + A . )!(1 _ A . }! ( 2 - 23 ) 

Ai=±l,0 A/v— il,0 i— 1 

Q °( 7 ) = 1 and Q^( 7 ) = for \n\ > N > (2.24) 

We note that the non-diagonal reduced helicity amplitudes T\.a±i = 1~\\±i are non-vanishing for any 
energy. As a result, the term ^ A (|'7a > a-i| 2 + |^a,a+i| 2 ) never vanishes, leading to a cross section that 
rises ~ at the threshold and contributing with a term ~ (1 + cos 2 6) to the angular distribution. Both 
elements differ clearly from the production of scalar particles. We therefore conclude that scalar spin-0 
particles in supersymmetric theories carrying muon-type charges can never be confused by fermionic 
charged particles. 

(b) Bosonic spectrum F[J = 1,2,...]: 

Restricting ourselves to CP-invariant theories, the electric and weak monopole charge term of any integer 
spin J tensor-field ip ai -- aj is accounted for by the current element 

3l = i K V ai - aj (2-25) 

It leads to P-wave production of the boson pair with the reduced helicity amplitude 

Tx\ = -^Qi(l) (2-26) 

Since the wave-function vanishes at the origin, the total production cross section rises ~ f3 3 at the 
threshold. Thus, opposite to wide-spread belief, the onset of the excitation curve near threshold does 
not discriminate the spin particle from higher integer spin J = 1,2,... particles. 

However, coupling the electroweak vector fields to the spin J fields in a consistent way [19], a non- 
vanishing magnetic dipole moment is generated for all particles with spin > 0. The non-zero magnetic 
current, which is proportional to 5 

j™ = -id»(^--- a '<p v]a2 ... aj ) (2.27) 

gives rise to non-vanishing off-diagonal helicity amplitudes which can be written, apart from an overall 
coefficient, as 

K x±1 = -^ ^^^J-i y^ [(^±A)Q^(7) + (JTA-l)Q^(7)] (2.28) 

The P-wave behavior near the threshold is reflected in the coefficient (3. The non-vanishing helicity 
amplitude T\ t \±i = T™^^ for J > is in apparent contrast to spin-0 scalars for which these amplitudes 



As before, asymptotic unitarity [19] modifies the coefficient of this current such that the gyromagnetic ratio is shifted 
again from g — J -1 to g = 2. 
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must vanish. Thus opposite to scalar production, higher spin J = 1,2,.. production will generate an 
additional term ~ (1 + cos 2 9) in the angular distribution, non-vanishing in the forward and backward 
direction. Thus the analysis of the angular distributions signals the zero-spin of the smuons unambigu- 
ously. 

Complementary to this theory-based argument on production properties, i.e. the onset of the ex- 
citation curve and the angular distribution, decay characteristics can also be exploited to supplement 
the analysis. The presence of the off-diagonal helicity amplitudes, Eq. (2.22) for fermions and Eq. (2.28) 
for bosons, implies that the final state particles F J and F J should be polarized when electron and/or 
positron beams are polarized. In addition, the relation Ai — A2 = ±1 for the non- vanishing off-diagonal 
helicities of the P£ pair produced via 7, Z exchange should lead to interesting polarization correla- 
tions which can be observed through the correlated decays of F,[ and F^ . In this case, the polar-angle 
distribution of the decay particles in the F J rest frame is described by the Wigner d function, 

V[F^a ai b a2 ] ~ d J Xa (9*) with a = a 1 -a 2 (2.29) 

where 9* denotes the polar angle between the F J flight direction and the ab axis in the F J rest frame. 
This configuration is realized in the dominant nx® decay channel of p,R. For scalar smuons c?q O (0*) does 
not depend on 9*. For J > 0, however, even if the final-state polarizations o"i,o"2 are summed over, the 
angular distribution is always non-trivial when the sum of two final-state daughter spins, j a and jb, is 
less than the spin J of the parent. In the opposite case, j a + jb > J, parity-violating decays in general 
guarantee non-trivial angular dependence 6 ; only in parity-preserving cases the decay distribution might 
be independent of the F J polarization. However, if the final-state polarizations a%, a 2 are measured, the 
9* dependence of d Xa (9*) for J > is always non-trivial, whatever values are taken for |A| and \a\ < J. 
After squaring the decay amplitude T>, the spin J can be determined by projecting out the maximum 
spin index 2J from the decay angular distributions. 

Therefore the analysis of the smuon decay distributions provides us with an alternative model- 
independent method for the determination of the zero smuon spin. 

2.4 Reconstructing the Event Axis 

(a) fi%jj<~jz in supersymmetry: 

The measurement of the cross section for smuon pair production fi%fi^ can be carried out by identifying 
acoplanar pairs [with respect to the beam axis] accompanied by large missing energy: 

e+e- - - (/i + x?)(^X?) - # (2-30) 

The analysis is model-independent and it provides unambiguously the onset of the excitation curve ~ /3 3 
near threshold. 

The construction of the production angle 9 is illustrated for the event topology in Fig. 3. For 
very high energy yfs S> mf 1R the flight direction of the daughter particles s can be approximated 

6 Only in exceptional cases, like r — > v T a\ with m T ~ s/2m ai , P-violating decays cannot be used as polarization analyzer. 
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by the flight direction of the parent particle [13] and the dilution due to the decay kinematics is small. 
However, at medium ILC energies the dilution increases, and the reconstruction of the fl R flight direction 
provides more accurate results on the angular distribution of the smuon pairs [12]. If all particle masses 
are known, the magnitude of the particle momenta is calculable and the relative orientation of the 
momentum vectors of fi^ and jl R is fixed by the two-body decay kinematics. The opening angles a± 
between the visible fi^ tracks and the parent £l r particles can be determined from the relation 

m 2 ± - m 2 ~Q = \fsE u ± (1 - 0~± cos a±) (2.31) 

The angles a± define two cones about the y + and yT axes which intersect in two lines - the true fx^ 
flight direction and a false direction. True and false solutions are mirrored on the plane spanned by the 
ji + and \X~ flight directions. Thus the flight direction can be reconstructed up to a 2-fold ambiguity. 




Figure 3: Event topology of the reaction e + e — > y\y R — ► M + Xi^ X? 

The characteristics of the angle 6 ft between the false and the true axis can easily be illustrated. 
If the decay planes of fl R and jl R coincide, the production axis is located in the common plane and 
the false axis coincides with the true axis. Rotating one of the two planes away from the other by an 
azimuthal angle </>, the angle 6 ft between the false and the true axis is related to <fi and to the boosts 
7± = 7(cosa^_ + [3)1 sina^_ with a*± being the ft R — > /U ± decay angle in the fl R rest frame with respect 
to the flight direction in the laboratory frame: 

2 sin 2 6 

cos 6 ft = l- -= ^ 2— (2.32) 

7+ + 71 + 27 + 7_ cos <p + sin (f> 

For high energies the maximum opening angle reduces effectively to 9ft < min(l/|7 + |, 1/|7_|) = 0(1/7) 

and approaches zero asymptotically when the two axes coincide. Quite generally, as a result of the 

Jacobian root singularity in the relation between cos8f t and <j>, the false solutions tend to accumulate 

slightly near the true axis for all energies. In total, the angular distribution of the false axis with respect 

to the true axis is given by 

dN V2 F [/3cos(0/i/2)] 



dcosOft y/l-cosOft 7 4 /? 4 (l-cos# /t + 2/7 2 /3 2 ) 2 



(2.33) 
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with F[0] = 1 at threshold and F[(3] ~ 7 for — > 1. The decrease of the coefficient ~ I/7 4 is 
compensated by the effective narrowing of 9ft in the denominator and by the increase of the function 
F for rising energy. Thus, the false axis is trailed by the true axis, mildly at low energies and tightly at 
high energies. Though the distribution of the false axis is flattened at low to medium energies compared 
with the original distribution of the true axis, the characteristic features are reflected qualitatively, 
nevertheless, cf. Fig. 4. For our theoretical investigation throughout the paper, the R-type slepton mass 
rri£ R = 300 GeV and sneutrino mass ra^ = 365 GeV are used and the chargino/neutralino mass spectra 
and the mixing elements are derived from the SUSY Lagrangian parameters Mi = 300 GeV, M\ = 150 



GeV, fj, = 500 GeV and tan/3 



10. This parameter set includes the lighter chargino mass m~± = 286 

Xi 



GeV and to the two lowest neutralino masses m^o = 148/286 GeV. 



1/2 
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Vj = 1 TeV 
m f = 300 GeV 
mj;=148GeV 
[arbitrary normalization] 
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UED : e e H fll H flI ^ H H Y,Y, 
reconstracted axis 
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Figure 4: (a) The angular distribution of the false reconstructed axis (blue shaded area) and of the ob- 
servable sum of true and false axes for smuon pair production in SUSY (black); (b) The same for 
pair production in UED. [The distribution of the false axis is slightly forward-backward asymmetric due 
to the parity violating fj,^ — > /i ± 7i decay.] 

Experimentally, the absolute orientation in space is operationally obtained by rotating the two /i^ 
vectors around the [V s1 axes against each other until they are aligned back-to-back in opposite directions. 
The flattened false-axis distribution can be subtracted on the basis of Monte Carlo simulations. For 
a fraction of events the production angle cannot be reconstructed, which in most cases is due to large 
initial state radiation and/or beamstrahlung reducing the nominal center of mass energy considerably. 
In addition losses occur due to measurement errors of the final particles. It is also important to note 
that background events, if they can be reconstructed at all under the wrong (mass) hypothesis, usually 
produce flat angular distributions and can thus be easily subtracted. 
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(b) fJ, R1 n R1 in UED : 

As proven in the previous subsection, the experimental observation of the sin 2 9 law determines the 
spin-0 character of new particles carrying non-electron fermion numbers, i.e. smuons, squarks, etc, un- 
ambiguously. This general conclusion can be illuminated by analyzing the spin-1/2 angular distribution 
of the UED KK excitation y^. The [non-normalized] angular distribution of Mj^Mri pairs in e + e~ 
collisions is described in Eq. (2.11): N(6) = 1 + cos 2 9 + (1 — 1 ) sin 2 6. The decay y R1 — > y ±r yi in the 
yf rest frame is governed by the right-handed coupling between the two leptons so that y~ , including 
the angle 9\ with the y R1 polarization vector, is preferentially emitted in the direction parallel to the 
y R1 polarization vector, D{0{) = (1 + kricos6>i) 2 with kr\ = ("i 2 lfll — 2m 2 1 )/(m 2 m + 2m 2 x ) ~ 1/3. 
The opposite rule applies to y R1 decays. 

Properly including the correlations among the two decay pairs, the predictions for the distributions 
of the true production axis, N(9), and the false axis, Nf(9), are displayed in Fig. 4(b). After subtracting 
the distribution of the false axis from the sum, the distribution of the true axis is markedly different 
from the distribution of the smuon polar production angle in Fig. 4(a). In particular, the production of 
spin-1/2 KK muons populates the forward and backward directions in contrast to spin-0 smuons. 

This method can be applied quite generally in e + e~ annihilation through 7 and Z exchange for any 
given theory. For the set {c\, C2} of coefficients in the true angular distribution 

N(9) ~ ci [1 + cos 2 6] + c 2 sin 2 9 (2.34) 

the false distribution Nf(6) can be generated unambiguously. Comparing the sum of the distributions 
of the experimentally reconstructed true and false events with N + Nj , the ratio of the two coefficients 
c<ijc\ can be fitted by using template methods. The fit will only be acceptable if at the same time the 
helicity structure of the decay vertex is chosen correctly. 

2.5 Experimental Analysis 

(a) Sparticle spectrum: In order to perform a realistic simulation of signal and background processes 
a sparticle spectrum is calculated using the program Isajet 7.74 [23]. The SUSY point described 
above can well be embedded in a mSUGRA scenario 7 with the parameters mo = 265 GeV,M]y 2 = 
375 GeV, tan /? = 10, Aq = 0, and sign/i = +1, corresponding to the Lagrangian parameters M\ = 
156 GeV,M2 = 291 GeV,/x = 488 GeV. The masses of the R/L-type smuons/selectrons, electron 
sneutrino, lighter chargino and two lightest neutralinos accessible at a 1 TeV ILC and relevant for the 
present experimental study are listed in Table 1. 

(b) Event generation: Events are generated with the program Pythia 6.3 [24] which includes initial and 
final state QED radiation as well as beamstrahlung [25]. The experimental simulation is based on the 
detector proposed in the Tesla tdr [26] and implemented in the Monte Carlo program Simdet 4.02 [27]. 

7 The particle masses corresponding to this reference point differ slightly at a level of a few GeV from the previously 
adopted masses in the theoretical illustrations. 
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Table 1: Spectrum of sleptons, charginos and neutralinos in the SUSY scenario (Mi = 156 GeV,M2 = 
291 GeV,/x = 488 GeVj. 
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The detector requirements are excellent momentum and energy resolution, good particle identification 
and full hermetic coverage. The detector response, resolution and particle reconstruction are treated 
in a parametric form. It is further assumed that the ILC can be operated at a flexible energy up to 
t/s = 1 TeV and that both lepton beams can be polarized at a degree of \V e - \ = 0.8 for electrons and 
\V e +\ = 0.6 for positrons. Beam polarization helps to enhance the production rates and to select the 
signal but it has no essential influence on the spin analyses of the distributions under investigation. 

(c) Event reconstruction : The reconstruction of the polar angle of pair production relies on the knowl- 
edge of the masses of the primary and secondary particles. Based on pure kinematics of two-body decays, 
like — ► i^Xii X\ ~ > W^Xi an d X2 ~^ ^x\i both masses can be determined from the energy spectra 
of the observable decay particle, see e.g. Ref. [28]. Alternatively, the excitation curve can be used to 
determine the mass of the primary SUSY particle pair. However, the observable cross section close to 
threshold is in general distorted considerably and the theoretical expectation has to be convoluted with 
initial state radiation (ISR), beamstrahlung (BS) and finite width effects. ISR can be rigorously treated 
in QED. The BS energy profile depends on the collider operation conditions, and it can be measured 
via Bhabha scattering [29]; it can also be calculated for given machine parameters [25]. The width of 
SUSY particles is calculable within a specific model, but can also be determined in a simultaneous fit 
of the excitation curve [16,17,30]. It can be safely assumed that the sparticle masses can be measured 
with a precision of one permille or better, see Ref. [31]. Such an accuracy is sufficient for the present 
study to reconstruct the event kinematics reliably. 

2.6 Simulation of e + e~ — > ArAr 

The detection of scalar smuons in the reaction e + e~ — > /i^/J^ with subsequent decays /x^ — ► n^Xi is 
relatively simple and clean. The energy spectrum of the decay muon is flat with minimal and maximal 
values given by E± = y/s/4 (1 — w| /m?) (1 ± (3) with (3 = (1 — 4m?/s) 1 / 2 . The event selection criteria 
are: (i) two oppositely charged ^ and nothing else in the detector; (ii) signed polar angle acceptance 
—0.90 < Q^cosOp < 0.75; where is the muon charge and the asymmetric cut rejects muons from 
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decays of W + W~ production; (iii) acoplanarity angle between the two muons A<J> W < 160°; (iv) the 
missing momentum vector should point inside the sensitive detector \oosQp | < 0.9; and (v) lepton 
energy within the kinematically allowed boundaries EL < En < E + (modulo resolutions). The resulting 
detection efficiency is typically around e ~ 0.60. 




Figure 5: (a) The unpolarized cross section of e + e~ — ► fJ^fi^ production close to threshold, including 
QED radiation, beamstrahlung and width effects; the statistical errors correspond to £ = 10fb _1 per 
point, the dash-dotted curve indicates, with the same coefficient, the hypothetical dependence a ~ (3 
instead of a ~ /3 3 ; (b) energy spectrum E^ from fi^ — > ^t~Xi decays; polar angle distribution cos9^ R 
(c) with and (d) without contribution of false solution. The simulation for the energy and polar angle 
distribution is based on polarized beams with (J> e -,V e +) = (+0.8, —0.6) at yfs = ITeV and C = 500 fb -1 . 
The smooth histograms represent high statistics expectations, the curve indicates a fit to the cross section 
(2.5). 

The unpolarized cross section as a function of energy close to the production threshold, including all 
instrumental effects, is shown in Fig. 5(a). The remaining flat background from W + W~ production (not 
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shown) amounts to about 0.5 fb, and it can be subtracted by extrapolation from the sideband below. 
The excitation curve exhibits a slow rise as expected from the characteristic dependence &fL R ji R ~ /? 3 
explained in Eq. (2.7). Such a behavior can be clearly distinguished from a much steeper hypothetical 
5-wave dependence a ~ /3, shown as well for comparison. 

The angular distribution is investigated in the continuum choosing an energy of yfs = 1 TeV and 
an integrated luminosity of L = 500 fb -1 . In order to enhance the signal and suppress SUSY and SM 
background processes, the beams are assumed to be polarized with values of (V e -,V e +) = (+0.8, —0.6), 
resulting in a cross section of 30 fb. 

The spectrum of the muon energy is shown in Fig. 5(b). The signal is very clean above the 
very low background from SUSY and W + W~ production, which gets reduced further to ~ 1.3 % after 
reconstruction of the kinematics. The primary flat energy spectrum, characteristic for a spin particle 
decay, is distorted due to acceptance cuts, event selection criteria and photon radiation. However, the 
minimal and maximal endpoint energies E± are clearly pronounced. The polar angle distribution of the 
smuons jlft, including both the correct and false solutions, is displayed in Fig. 5(c). The tiny background 
gives a flat contribution. The false solution exhibits a large pedestal with some enhancement in the 
central region reflecting mildly the primary distribution (cf. (5 = 0.8). The ambiguous solution can 
be calculated using Monte Carlo simulation and it is subtracted in Fig. 5(d). The expected sin 2 9 
distribution of Eq. (2.8) is clearly visible. A fit of the shape to the experimental angular distribution 
yields 

cfcr exp 

~ 1 + a cos + 6 cos 2 6> (2.35) 

a cos 

a = -0.020 + 0.016 and 6 = -0.979 ± 0.022 

confirming the conjectured forward-backward symmetric sin 2 behavior of spin-0 smuon production 
with high precision. 



3 SELECTRONS 

3.1 Production Channels in e + e~ Collisions 

In the selectron production process the lepton number can flow from the initial to the final state. 
Therefore, besides the e + e~ annihilation channel mediated by 7, Z exchange, cf. Fig. 1(a), also t-channel 
exchange of neutralinos, cf. Fig. 1(b), contributes to some channels. [In e~e~ collisions selectrons are 
produced solely by t-channel and u-channel exchanges.] Among all these channels, which are summarized 
comprehensively in Table 1 of Ref. [17], pair production of e^e^ is easiest to control experimentally 
if, as realized in many models, the R-type selectron is significantly lighter than the companion L-type 
selectron. Equal-particle channels are also preferred theoretically; their analysis is most transparent by 
including the standard annihilation which is well controlled. 
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Two electron/positron polarization states can generate the e~^e R pair: 

e R e L ~> ZrZr bi Z exchange] (3.1) 

e l e R -> CrCr b, z >X° exchange] (3.2) 

Though the signal process (3.1) is the analogue of smuon pair production in e + e~ annihilation, we cannot 
anticipate that in rival J = 1/2, 1, .. processes t-channel exchanges do not occur when the lepton number 
can flow from the initial to the final state. Moreover, even electron polarization [32] can be realized 
only at a degree < 1 so that impurities from the process are mixed in in any case. It is thus plausible to 
evaluate e R e R pair production for unpolarized electron/positron beams. [The general expression of the 
polarized cross section is given in the Appendix.] This case exemplifies all the interesting characteristics. 

After exploiting the conservation of the lepton current, the spinorial parts of the matrix elements 
for s-channel 7, Z exchange and f-channel neutralino exchange are identical, ~ {v e +l^, u e-) ^ with fc M 
denoting the (spacelike) 4-momentum transfer. The t-channel contribution can therefore be mapped 
into generalized charges, introduced in analogy to Eqs. (2.2) and (2.3): 

Ql = 1 + {s 2 w - 1/2) c^D z (s) (3.3) 

4 

Qr = 1 + t 2 w D z (s) + iNkifc^D^it) (3.4) 
fc=i 

The indices R and L denote the electron helicities. The pole part of the xt propagator, k = 1, ..,4 in 
the Minimal Supersymmetric Standard Model (MSSM), is denoted in the center-of-mass frame by 

°<-." |5 rk ^>-^« with ^ = 1 - 2 «-'4;)/ s < 3 - 5 > 

while N denotes the neutralino mixing matrix, see Ref. [33]. Near the threshold, the e^e R pair is 
produced in a P-wave with amplitude ~ (3. With rising energy however an increasing number of orbital 
angular momenta is excited and the propagator starts diverging in the forward direction ~ s/m 2 [for 
rrie R > m x° k after running through a maximum at cos# ~ Afc//3g H ]. 

The differential and total cross sections can be cast into the form 

d A_ [e + e -->g£eg = sin 2 [ Q \ + Q \\ (3.6) 

" 2 /3 3 [<Q|sin 2 0) + (Q|sm 2 0}] (3.7) 

with {Q 2 L R sin 2 0) = f dcos 9 sin 2 Q 2 L R . Mass and energy dependence of the integrated charges can be 
adopted from Ref. [17]: 

<Qisin 2 0> = ^ [1 + {s 2 w - 1/2)6^ D z {s)} 2 (3.8) 

. 4 

(Q 2 R sm 2 9) = -[l + t 2 w D z (s)] 2 + 8[l+t 2 w D z (s)}c^Y,\ N n\ 2 f j 

4 3=1 
+8c^ \NjiN kl \V (3.9) 
j,k=i 
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with the coefficients 



A 2 -B 2 \- + R [ -2/3 + Aj-ln^l j = k 

f = -PA J + %-^ln^±l and = \ ' ~ (3.10) 



It follows that the production of e-type super symmetric scalar particles is characterized by the 
following two rules: 



#1 threshold excitation 
#2 angular distribution 



P 3 

sin 2 (9 e?(cos<9) 

sin 2 9 near threshold 



(3.11) 



(3.12) 



Independent of energy, the angular distribution must behave ~ sin 2 9 close to the forward and backward 
directions where it must vanish by angular momentum conservation. While this behavior may be 
masked in practice by the singularity in Q developing in the forward direction at high energies, no such 
interference will arise in the backward direction. Since the Xfc exchanges give rise to a P-wave near 
the threshold, in the same way as 7, Z exchange, a simple picture with a ~ /3 3 and do jd cos 9 ~ sin 2 9 
emerges at the threshold. 
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Figure 6: (a) The threshold excitation for selectrons; and (b) the angular distribution in e + e~ — > e^e^, 
for the SUSY parameters specified in the text. These results are compared with the production e + e~ — > 
e ~Ri e Ri °f elect™™ Kaluza-Klein states in UED. The first KK mass Mk 1 is taken to be 295 GeV. 



Asymptotically the total cross section scales as 

ira 2 1 s 



a 



for 



:3C 



c w s 



(3.13) 



as expected from the forward enhancement of the i-channel exchange. 

These characteristics are displayed quantitatively in Fig. 6(a/b) for the onset of the excitation curve 
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and the angular distribution close to threshold. 



3.2 KK Excited States ef in UED 

The analysis presented above repeats itself rather closely for the KK excited states carrying electron 
lepton number; again we choose the first KK excited R-type electrons with vector couplings to not 
only 7 and but also to Z bosons as a representative example: 



e e 



e m e m 



(3.14) 



Analogously to Fig. 1(b), the t-channel exchange of the vector and scalar KK excitations B\ and b\ [the 
supplement left with B\ from the 5-dimensional vector] add to the standard 7, Z exchanges corresponding 
to Fig. 1(a). 

Despite the complicated superposition of vector and scalar interactions, Fierzing techniques allow 
us to cast all contributions into the s-channel 7^ (8) 7^ form 8 for the chiral a, j3 = L, R elements: 



M[e + e -> e^e ffl ] = Q l af} [v(e + ) ltl P a u{e )] [u{e Rl )^Ppv{e + Rl )\ 

with the bilinear charges [34]: 

Qll = 1+ {sw - 1/2) c^D z (s), Qi R = Ql L - l J ^- T ? w D bl (t) 

Qrl = 1 + t 2 w D z (s) + c ^^D Bl (t), Q RR = 1 + t 2 w D z (s) + Cw 2 D Bl (t) 



(3.15) 



(3.16) 



Apart from standard notations, denotes the first KK mass, and the t-channel propagators are 

defined as D Bl (t) = D bl (t) = s/[t - M 2 Ki ] = -2/[l - 2(m 2 eRi - M 2 Ki )/s - em cos0]. 

After introducing the familiar quartic charges 



^1 = 4 [\Qrr\ 2 + \Qll\ 2 + \Qrl\ 2 + \Qlr 2 
Ql = % [QrrQrl + QIlQlr] 



Q\ = l [\Qrr\ 2 + \Q\l\ 2 - \Qrl\ 2 - \QW 2 ' 

the differential cross section can be written in the compact form 



da 



net 



e e 



e Rl e Rl\ 



(1 + P 2 cos 2 9)Q\ + (1 - p 2 )Ql + 2(3 cos 9Q\ 



dcosO L ~ ' 2s 
from which the total cross section follows by integration over the polar angle: 

„2 



cr e e 



e Rl e RV 



ira 
~2s 







(Q\) + p 2 (cos 2 9Q\) + (1 - (3 2 ){Q\) + 2/3(cos 9Q\) 



(3.17) 



(3.18) 



(3.19) 



The matrix element e + e — > /J-R 1 fJ, R1 has earlier been denned implicitly in the same way, the charges Qr and Ql to 
be identified with Qr = Qrl = Qrr and Ql = Qll = Qlr in this case. 



17 



Both observables can serve as discriminants for Kaluza-Klein states against super symmetric selectrons. 

By inspecting the cross sections in Eqs. (3.19) and (3.18) we can easily conclude, without studying 
details, that 

a[e + e~ — > e^e^-J ~ j3 near threshold (3.20) 

[e + e-^e+ iem ] - (1 + [3 2 cos 2 9) £(cos 6) + • • • 

— ► flat in cos 9 near threshold (3.21) 



dcos 9 



as generally expected for fermion pair production near the threshold. As for smuon pairs, these results 
contrast strongly to supersymmetric scalar e production. Most striking is the non-vanishing angular 
distribution in the forward and backward directions. This is exemplified quantitatively in the comparison 
of Fig.6(a/b). 

Asymptotically however the total cross section, unlike the previous examples, approaches a non-zero 
value 

ttq 2 1 , 

for s -> oo (3.22) 



r 4 M 2 
C W 1V1 K 1 

due to the enhancement in the forward direction, which is a remnant of the Rutherford pole damped 
by the Yukawa mass cut-off in the exchange of heavy particles. 



3.3 General Analysis 

Independent of the lepton number flow by additional i-channel exchange mechanisms, the s-channel 
7, Z exchange in the production of charged fermions of any spin J = 1/2, 3/2, . . . will generate the pair 
in an S'-wave so that the cross section should rise at threshold ~ @ in contrast to the scalar /3 3 particle 
production. The same 7, Z exchange mechanism will generate a non- vanishing angular distribution in 
the forward and backward directions 7^ sin 2 9. 

In contrast to the production of muon-type pairs, the additional i-channel exchanges in the pro- 
duction of J = 1, 2, .. integer spin electron-type pairs will in general give rise to an 5-wave component 
~ (3 in the onset of the excitation curve. Since all spin J > particles in asymptotically well-behaved 
field theories [19] will carry a non-vanishing magnetic dipole moment, the angular distribution for both 
muon-type and electron- type pairs will not vanish in forward/backward direction. This argument can 
be supplemented by studying the polarizations ~ d^ff) in the decays of the two spin J particles. 

Thus in parallel to the smuon case, also for selectrons in supersymmetric theories experimental paths 
can be designed for establishing the scalar spin-0 character unambiguously. 



3.4 Simulation of e+e — > ^r^r 

The detection of scalar selectrons in the reaction e + e~ — ► e^e^, — ► e + Xi e ~Xi is again very clean. The 
event simulation, selection and analysis proceed in complete analogy to smuon production, described in 
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the previous section, just replacing the observable leptons by an e + e pair. 
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Figure 7: (a) The unpolarized cross section of e + e~ — > e^e^ production close to threshold, including 
QED radiation, beamstrahlung and width effects; the statistical errors correspond to £ = 10 fb" 1 per 
point, the dash-dotted curve indicates a hypothetical dependence a ~ j3; (b) energy spectrum E e from 
— > e ± x? decays; polar angle distribution cos6e R (c) with and (d) without contribution of false solution. 
The simulation for the energy and polar angle distributions is based on polarized beams with (V e - , V e + ) = 
(+0.8,-0.6) at -^/s = 680 GeV and C = 200 fb -1 . The smooth histograms represent high statistics 
expectations, the curve indicates a fit to the cross section (3.6). 

The unpolarized cross section at threshold is displayed in Fig. 7(a). The expected event rates are 
larger than for smuon production due to additional i-channel neutralino exchange. One observes a cross 
section typical for P-wave production of spin-0 particles with a dependence a-+— ~ /3 3 , as explained 

e R e R 

in Eq. (3.11). The excitation curve may be easily distinguished from a much faster rising hypothetical 
a ~ (3 behavior, shown for comparison as well. 



The study of the e~J,e R production is performed close to threshold in order to separate out, as well 
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as possible, the factor sin 2 9 in the polar-angle distribution, see Eq. (3.12). An energy of y/s = 680 GeV 
appears to be a good compromise between signal and background event rates. An enhanced signal is 
obtained by choosing beam polarizations of (V e -,V e +) = (+0.8,-0.6). The cross section amounts to 
47 fb, and an integrated luminosity of £ = 200 fb" 1 is assumed for the simulation. The expected electron 
energy spectrum E e of the decays — > e^x? is shown in Fig. 7(b). It exhibits a clear signature above 
a negligible background from VF + VF~ and e^ej production. 

The angular distribution cos 9i R , shown in Fig. 7(c), is still peaked towards the forward direction, 
due to the remnant i-channel x° contributions, above a fairly constant pedestal from the false solutions. 
The spectrum after subtracting the ambiguity, displayed in Fig. 7(d), vanishes in the very forward and 
backward directions, reflecting the overall sin 2 9 factor. A fit according to the differential cross section 
formula (3.6) yields a very good description of the simulated data. As a by-product, the results of the fit 
can be used to determine or cross-check the neutralino mixing parameters | iVfc 1 1 2 entering the expression 
(3.4) of the generalized charge Qr. 




cos0„~ cos0 P . 



Figure 8: Polar angle distribution cos#g fl of e\e~ L — ► (a) with contributions from background 

and false solution; (b) after subtraction of background and false solution. The simulation is based on 
polarized beams with (V e -,V e +) = (—0.8, +0.6) at \/s = 680 GeV and £ = 300 fb -1 . The smooth 
histograms represent high statistics expectations, the curve indicates a fit to the cross section (3.6) 

The i-channel neutralino exchange can be considerably reduced by choosing opposite beam polar- 
izations of (V e - ,V e +) = (—0.8, +0.6). These conditions however imply a much lower cross section of 
a = 3.4 fb at the same energy y/s = 680 GeV and a significantly larger background from WW and e~Re~L 
production. The results of simulations assuming £ = 300 fb^ 1 are displayed in Fig. 8. The polar angle 
distribution is much flatter and shifted towards the central region, approaching the expected sin 2 9 law 
for completely polarized beams. A fit to the subtracted spectrum exhibits a skewed sin 2 9 distribution, 
reminiscent of small, residual t-channel contributions. 
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4 CHARGINOS AND NEUTRALINOS 



4.1 Production Channels in e+e Collisions 

The prototypes of non-colored super symmetric spin-1/2 fermions are the charginos xt an d the neutrali- 
nos Xi an d X2- They are produced in diagonal and mixed pairs in e + e~ annihilation: 

e + e~^xtXi with xf^ff'x ! (4.1) 

e + e~ - xixl and xlxl with xl - fM (4-2) 

Though a significant fraction of the decays is mediated potentially by f intermediate states as predicted 
in the reference scenarios SPSla/la' [18,35,36], other decay modes can still play a significant role due 
to large production cross sections, in particular for diagonal pairs. The mixed neutralino production 
channel X1X2 is easier to analyze in the threshold region when studying the onset of the excitation curve, 
but the diagonal pair X2X2 gives rise to a better textured visible final state that allows the reconstruction 
of the flight axis up to a 2-fold ambiguity, while the axis can be reconstructed for mixed pairs only if 
X2 cascades down through an intermediate slepton [37]. 




Figure 9: Two mechanisms contributing to the production of X1X1 an d three to X1X2 an( ^ X2X2/ 
s-channel "y/Z exchanges and t/u-channel v and e exchanges, respectively. 



Two mechanisms contribute to the production of xfxi anci three to X1X2 an< ^ X2X2- s-channel j/Z 
exchanges and t/u-channel v and e exchanges, respectively, as illustrated in Fig. 9. 

By using Fierzing techniques, both the i/u-channel diagrams can be mapped onto the s-channel 
diagram, generating the bilinear charges Q a p [33,38]: 



X1X1 



Ai Aj 



Qll = D l -F l cos 2<\> l Q LR = D' L -F' L cos 2<t> R 
Qrl = D R - F R cos 2(p L Q RR = D R - F R cos 2cj) R 

Q LL = +AZ ijS -^{s 2 w - l/2)D z {s) - g Lij D iL (u) 
Q LR = -4Z* jS -^(s* w - l/2)D z (s) + g* Lij D, L (t) 
Qrl = +2ijC^D z (s) + Q Ri jDe R {t) 
QRR = -2* jC ^D z (s) - g* Rij De R {u) 



(4.3) 



(4-4) 
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[in the usual notation s?, w = sin 2 29w]. The normalized propagators in these charges read 

D L = 1 + As^(s 2 w - l/2){s 2 w - 3/4)D z (s) F L = 4s^(s 2 w - l/2)D z (s)/A 
D R = 1 + c^(s 2 w - 3/4)D z (s) Fr = c w 2 D z (s)/A 

D' L = D L + s^D 9 (t)/A F' L = F L - s^D,(t)/A (4.5) 

and <f>L,R denote the mixing angles rotating the gaugino/higgsino current to the chargino mass basis; the 
rotation angles are determined by the SUSY Lagrangian parameters M 2 , ^ and tan (3 [38]. The matrices 
Z and Ql,r are combinations of the mixing matrix elements N in the neutralino sector [33]: 

Zij = (N i3 N* 3 - N i4 N* 4 )/2 

QiAj = {N i2 cw + Nnsw)(Nj 2 cw + Nj 1 s w )/s% w 

Qmj = NaN^/cyy (4.6) 

They are derived from the Lagrangian parameters noted above and supplemented by the U(l) gaugino 
parameter M 1 [in the MSSM] . 

Defining the quartic charges Qi,Q 2 , Q3 in the same way as Eq. (3.17), the differential and total cross 
sections can be written as 

* (3{[1- ((4 - ii 2 -) 2 + /3 2 cos 2 0]Qi + 4 W Q 2 + 2/JCOS0Q3} (4.7) 



dcosO 2s 

na 2 f. 



a 



■P{[l-( f i 2 - f j 2 )](Q 1 )+f3 2 (cos 2 6Q 1 )+4 t i lf i j (Q 2 ) + 2(3(cos8Q 3 )} (4.8) 

generically for any pair of masses with /ij = rrii/y^s, and f3 2 = [1 — (fa — //j) 2 ][l — (fa + fJ-j) 2 ] coinciding 
with the velocity squared for equal masses; f s = 1 or 1/2 denotes the statistics factor for pairs of unequal 
and equal particles, respectively, in the final state. 



(a) Charginos: 

Near the threshold the cross section rises ~ (3 since the charged Dirac particles are generated in S'-waves. 
For asymptotic energies the cross section scales as 

47TQ 2 

a — — (1 + Aj±) for s -> 00 (4.9) 

oS 

complemented by the coefficient A^± = 0(1) including mixing matrix elements. 

The angular distribution near the threshold is flat and need not vanish as for scalar particles in the 
forward/backward direction. With rising energy the sneutrino t-channel exchange excites an increasing 
number of higher orbital angular momenta and thus modifies the familiar spin- 1/2 asymptotic distribu- 
tion to ~ (1 + cos 2 9) G (cos 9). 

The characteristics for supersymmetric chargino production can therefore be summarized in the 
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following points: 



#1 threshold excitation ~ j3 (4-10) 
#2 angular distribution ~ (1 + j3 2 cos 2 9) £?(cos 9) + ■ ■ ■ 

— > isotropic near threshold (4-11) 

As will be argued later, these two characteristics can be mimicked by higher half-integer spin states. 
Thus, the observation of the characteristics (4.10) and (4.11) is necessary for chargino spin assignments 
in supersymmetric theories but not sufficient. The production characteristics must be complemented 
by decay characteristics to determine the spin of charginos unambiguously. 



(b) Neutralinos : 

The Majorana nature of the neutralinos forbids the 5"-wave production of the diagonal X2X2 P arr at 
threshold with equal spin components along the e + e~ beam axis as a consequence of the Pauli principle. 
This conclusion can also formally be drawn by observing that near threshold the sum of the quartic 
charges is reduced to Qi + Q2 = [\Qll + Qlr\ 2 + \Qrl + Qrr\ 2 ]/^ so that the final-state current 
becomes purely vectorial, forbidden however for neutralino Majorana fields which can only be coupled 
to axial- vector currents. The P-wave production mode leads to the onset of the excitation curve ~ /3 3 . 
The angular distribution follows the spin-1 rule for 7/Z exchange [39], modified however by a spin-1 
and spin-0 mixture from selectron t/u exchanges. Inserting the quartic charges in Eq. (4.7), the angular 
distribution is given near the threshold by 

dathr 



d cos 9 As 

k=L,R 



Z£(l + cos' 9) + Ag k (2g k - Z k ) sin J 9 + %5 k g k {5 k ~g k + Z k ) cos 2 9 



(4.12) 



where the coefficients are defined in terms of the matrices Z and Ql,r as 

Zr = Z<$D z 9r = Gr[1 + 6 r ] (4.13) 

Z L = 4Zs~l(l/2 - s 2 w )D z g L = G L [l + S L ] (4.14) 

with 6l r = (mlo — m| )/(m 2 + vn\ ). For large selectron masses the s-channel Z contributions are 

' X2 L,R X2 L.R 

dominant and the angular distribution is reduced to (1 + cos 2 #), characteristic for Majorana fermion 
pair production. However, if for particle masses of the same size, the t— and u— channel selectron 
contributions are dominant, the angular distribution is in general a mixture of sin 2 9 and cos 2 9 terms 
with coefficients varying with the particle masses. Above the threshold, higher orbital angular momenta 
are excited by the selectron exchange mechanisms, not altering however the asymptotic behavior 



47ra 2 
"3s" 



<7 — > — - — A^o for s — > 00 (4-15) 



with the coefficient A^o = 0(1). 

Also mixed X1X2 P a hs will be produced near the threshold in a P-wave if their CP parities, H, are 
equal. If they are different however 5-wave production is possible and the cross section rises ~ 0. In 
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theories with CP violation S- wave production is predicted in general [33,40]. 
These observations are summarized in the following rules: 

$d threshold excitation ~ /3 3 for X2X2 (4.16) 

~ /3 3 //3 for X1X2 [ident./diff. Major ana phases] 
#2 angular distribution near threshold ~ cos 2 / sin 2 9 mix (4-17) 

These points are illustrated for charginos and neutralinos in Figs. 10 and 11, respectively. The 
parameter set introduced earlier, gives rise to the chargino mass m-± = 286 GeV and the neutralino 
masses m x° /2 = 148/286 GeV. The residual linear /3-dependence, = 0.14, generates a slight increase of 
the angular distribution with cos 6. For neutralinos the chosen parameter set leads to a dominant sin 2 9 
component in the angular distribution, supplemented however by small additional contributions, cf. 
Eq. (4.12), which render the distributions non- vanishing at the very edges of the forward and backward 
directions. These results will be confronted with phenomena in UED and a general analysis in the next 
subsections. 
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Figure 10: (a) The threshold excitation for charginos and (b) the angular distribution in the process 
e + e~ — > XiXi f or the SUSY parameters specified in the text; both compared with W^W^ pair production 
in UED. 



4.2 KK Excited States W* and Z 1 /j 1 in UED 

The counterpart of X\ an d X2 * n UED are the KK excitations and Wf = Z\ , while x± an d B\ = 71 
are the stable particles of the two theories with minimal mass 9 to which all other particles cascade down. 

9 Electroweak mixing at the KK level is neglected as before. 
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Figure 11: (a) The threshold excitation for neutralinos and (b) the angular distribution in the processes 
e + e~ — > X1X2 an d X2X2 f or the SUSY parameters specified in the text; both compared with Z\Z\ pair 
production in UED. 
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Figure 12: Two mechanisms contributing to the production of W 1 W 1 and two to Z\Z\; s-channel 
7/Z exchanges, t-channel v\ and t/u-channel en exchanges, respectively. 

The cross sections for the processes 

e+e" -» WfW{- (4.18) 
e+e" -» Z X Z X (4.19) 

are closely related to the corresponding SM processes e + e~ — ► W + W~ and ZZ, cf. Fig. 12. Standard 
SM couplings are attached to the currents, and the exchanged neutrino and electron must be substituted 
by the heavy KK excitations. In the limit in which the masses of the t/w-exchange leptons are neglected, 
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the cross sections approach the SM form of Refs. [41,42]. 

The differential and total cross sections for can be expressed by the generalized charges 



Ql = 1-(s 2 w - l/2)s^D z {s) and Q R = 1 - D z (s) 



In this notation they can be written as 
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dcos 9 
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~8s~ 



P 
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s^F x {s,9) + -(Ql + Ql)F 2 ( S ,6) - Q L s^F 3 (s,9) 
s w + \(Ql + Or) ^(s) - Qlsw °3 (s) 



The angular functions F(s,9) and the energy-dependent coefficients a(s) are given by 

7 2 (1 + p 2 - 2(3 cos 9) 2 (4 + 7 2 /? 2 sin 2 6) + 4(3 2 sin 2 9 
>~ 4(A-(3cos9) 2 

F 2 {s,6) = 2(3 2 [4 7 2 (4 + 7 2 /? 2 sin 2 #) + 3sin 2 #] 

2 7 2 /? [2(3 - (1 + (3 2 ) cos 6} (4 + 7 2 /? 2 sin 2 0) + 3/? 2 sin 2 



(A - (3 cos i 



(4.20) 

(4.21) 
(4.22) 

(4.23) 
(4.24) 
(4.25) 



P 2 ) 



with 7 = y/s/2m w ±, A = 1 - 2{m 2 w ± - mlj/s, and 

ai (a) = 4( 7 4 + ll 7 2 - 3)/3 -(2 7 2 -l)r ? 2 -2 ?? 4 -2[7- 2 + 2(2- 7 - 2 )7 ? 2 ]/(A 2 

+ [4A - (7 + A)r/ 2 + (1 + 2A)t? 4 ] L{(3)/A 
a 2 (s) = 8/? 2 (4 7 4 + 20 7 2 + 3)/3 

a 3 (s) = 6A + 8(2 7 4 + 9 7 2 - 5)/3 - [ 7 ~ 2 + (2 - 7~ 2 )? ? 2 ] (2 7 2 - 1 + r/ 2 ) 

- {6(A 2 - (3 2 ) + [ 7 " 2 + (2 - 7" V] [7 + A - (1 + A)?? 2 ] } L(fi)/A 

with 77 = m^/m^yi, (3 = (3/ A and £(x) = (l/2x)log[(l + x)/(l — x)] — > 1 as x — > 0. Although each of 
the individual coefficients cti(s) grows as s, unitarity cancellations reduce the sum of all contributions 
to the expected scaling behavior of the cross section [41,42]: 
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The logarithmic term is generated by the KK neutrino exchange mechanism. 
For Z\Z\ production the differential and total cross sections read 
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with r] = nie^/mzj, A = 1 — 2{m 2 z ^ — m 2 Li )/s and 7 = 7 A. For asymptotically large energies, the 
standard behavior 
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is predicted for the total cross section. 

Near the thresholds the total cross sections rise as 



a[e + e~ -> W^W^ , Z ± Z ± ] ~ (4.33) 



while the angular distributions 



-Al^ll - l + 0(f3)cos8 and - -^—[Z 1 Z 1 ] ~ I + 0(/3)cos 2 # (4.34) 
<t a cosp 2 a a cost/ 2 

are essentially flat in the threshold region. The flat behavior is modified however linearly in (3 above 
the threshold as evident from Fig. 10(b). 

Comparing the predictions for the spin-1 KK excitations of the weak gauge bosons with the spin- 1/2 
charginos and neutralinos, we arrive at a mixed picture, cf. Figs. 10 and 11. In the chargino sector the 
onset of the excitation curves does not discriminate one from the other. However, due to the Major ana 
nature of the neutralinos, the onset for X2X2 ~ ^ ls different from Z\Z\ ~ (3. 

Final state analyses are necessary to discriminate charginos from KK W 1 bosons. Due to the 
vectorial/axial- vectorial couplings, in both theories, the electron-positron pair annihilates in a spin-1 
state polarized parallel to the beam axis. Angular momentum conservation then demands the same 
polarization state for the charginos which are coupled in an S-wave. Choosing longitudinally polarized 
electron beams with a degree close to one [32], the decay angular distribution is dictated by 

mi ~+ (//')X?] ~ 4 7 >*) ~ cos(0*/2) or sin(0*/2) (4.35) 

[depending on whether the initial x± helicity A and the difference a = cr(xi) ~ a {ff) °f final-state 
helicities are of equal or opposite sign]. This can be contrasted to the polarization of the W-f" which 
must be either 1 or 0, so that, in the same notation as before, 

V[Wi -»■ (/f)7i] ~ d\ a {6*) ~ (licosfl*) or cos6>* or sin^* (4.36) 

[depending on whether |A| = \a\ = 1,0 or otherwise] with quite a different Wigner d function compared 
to the supersymmetric signal. Thus the final state analysis provides a clear discrimination. 

In conclusion. The supersymmetric chargino/neutralino sector can be discriminated in spin analyses 
from the KK excited weak-boson sector in theories of universal extra space dimensions, but final state 
analyses of the decaying particles are required. 



4.3 General Analysis 

(a) Charginos: S'-wave production of chargino pairs gives rise to the (5 onset of the excitation curve near 
the threshold. This behavior is expected for all charged half-integer spin Dirac particles. In parallel, the 
angular distribution in the production process does not discriminate the particles. Bosons with spin > 1 
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also follow the S-wave pattern if they are produced pairwise through t- and/or u-channel exchanges. 

Quite generally, the (polarized) electron/positron pair either annihilates in a polarized spin-1 state 
for vector currents, as exemplified above. 10 This gives rise to polarization effects d J x {9*) in the F J decays 
and to correlations of the form 

d{±i,M)d{A0*2) and/or d J X(7 (9l) d{ a ,(9* 2 ) (4-37) 

between the angular distributions of the decay products of the particle pair F J F J which is generated 
in an 5-wave near threshold. The characteristic dependence of the d functions on J can be exploited to 
determine the spin. [Details will be presented in the subsequent experimental subsection.] 

(b) Neutralinos : For clarity we focus on the production process (2.13) for equal-type particle-antiparticle 
pj .J. pJ anc j p ar ticle-particle F = F J pairs. As argued before, S"-wave production is expected in 
general if the neutral fermions F J and F J are different from each other; it gives rise to the (3 dependence 
of the cross section near threshold as opposed to the /3 3 production law of the Majorana particle x°>. 

It has been shown quite generally in Ref. [39] that Majorana pairs F J F J are always produced in P 
waves near threshold, with a (1 + cos 2 9) angular distribution for spin-1 j/Z exchange. If t/tt-channel 
exchanges are switched on, S"-wave production of Majorana fermion pairs remains suppressed for all 
interactions conserving electron-chirality. While the rise of the excitation curve ~ /3 3 does not change, 
the angular distribution is modified however to a mix of cos 2 6 and sin 2 9 terms. 

Thus the spin of charginos and neutralinos cannot be discriminated unambiguously unless the stan- 
dard correlation tests involving the chargino/neutralino decays with reasonable polarization analysis 
power are performed. The analysis of polarization effects in F J decays, eventually supplemented by 
correlation effects in double F J F J and F J F J decays, in the way discussed above, will lead to the un- 
ambiguous spin assignment J = 1/2 of the xf chargino and the X2 neutralino. 

4.4 Simulation of e + e~ — > xtx.i an d e + e~ — ► X2X2 

Chargino production and detection proceeds via e + e~ — > x±Xi ~~ * W + X? W~x\ with a branching 
ratio B(x\ ~> W^Xi) = 1 in the reference point considered. Distinct experimental signatures are 
either purely hadronic decays WW — > qq' qq' — > 4 jets or mixed hadronic and leptonic decays WW — > 
qq' Iv — > 2 jets + 1 lepton. For a complete reconstruction of the kinematics, including production and 
decay angles, only the 4-jet final state can be used. However, information on the individual W^ charge 
is heavily spoiled by large fluctuations during the fragmentation process which may lead to track losses 
and/or wrong track assignments to the parent particle. Only a folded angle |cos#| can be obtained. 
In contrast, the electric charge of individual W^ can be identified in the mixed hadronic and leptonic 
decays, W±W^ qq'i^v. 

A potential background source is neutralino production e + e~ — ► X2X2 with subsequent decays X2 ~^ 

10 Scalar and tensor couplings, which would correspond to spin-0 states, vanish in the limit of zero electron mass due to 
electric chirality conservation. 
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Zxi (B = 0.13) and x\ ~ > ^X? (B = 0.87). The hadronic decays of Z and h provide an event topology 
and kinematics very similar to chargino production. A distinction may be possible on the basis of 
excellent di-jet mass resolution as anticipated in the design of future ILC detectors [43]. The goal is to 
achieve an efficient separation of hadronic W and Z decays, which also implies a reliable identification 
of the heavier Higgs decays. 




Figure 13: Cross sections of (a) e + e~ — > XiXi production and (b) e + e~ — > X2X2 production close to 
threshold, including QED radiation, beamstrahlung and width effects. The statistical errors correspond to 
C = 10fb _1 per point, The dash-dotted curves indicate hypothetical cross section dependencies ~ 
f3 3 and c^,o^o ~ j3 for charginos and neutralinos, respectively. 

The cross sections for X\X\ an d X2X2 production as a function of energy are shown in Fig. 13(a) 
and Fig. 13(b), respectively. Since the masses are almost degenerate, the threshold energies of both 
reactions are very close and practically coincide. However, the chargino cross section rises much faster 
with a ~ [3 compared with the slow onset of the neutralino excitation curve a ~ /3 3 . It is obvious from 
the threshold curves that the different (opposite) (3 dependence for the two reactions can be easily ruled 
out. 

As pointed out in the previous section, polarization effects in xf and X2 decays must be exploited 
to determine the spin J = 1/2 of the xf chargino and x\ neutralino unambiguously. 

(a) Charginos xf'- 

The charginos in the production process e + e~ — > X.1X1 are polarized and even the polarization averaged 
over the production angle 9 is in general non-zero (also for unpolarized beams). The cosine of the decay 
angle 0\ v ± between the momentum direction in the chargino rest frame and the X\ momentum 
direction in the laboratory frame, identical with the spin quantization axis, can be determined by 
measuring the VF ± energy in the hadronic decay — > qq': 

E w± = 7 (E* w± + 13 p* w± cos 9* w± ) (4.38) 
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where 7 = ^/s/2m~± and j3 = (1 — 4771^/s) 1 / 2 . The Vl^-boson energy and momentum in the chargino 
rest frame, E^± and p^±, can be derived from the xfiXi an d masses. Furthermore, the electric 
charge of the individual and xt can be identified by tagging the electric charge of the from 
the other chargino x^f decay through the leptonic mode — > f^vi. All these features can be used to 
study the x^-polarization through the angular distribution in the two-body decays Xi ~~ > W^X?- For 
the chargino as a spin- 1/2 particle the decay distribution is linear in cos^-t: 

where the coefficient («jy±) is the product of the %i polarization averaged over the production distri- 
bution and the polarization analysis power of the decay mode Xi ~~ * W^Xii f° r details see Ref. [44]. 
The two coefficients (k w ±) are identical as a consequence of CP symmetry. The average helicities of X\ 
and Xi have the same magnitude but opposite sign in e + e~ annihilation for arbitrary beam polarization. 
For non-zero (k^±) the decay angular distributions provide a unique signal for the spin J = 1/2 of the 
chargino xf . 

The decay angular distribution is studied in an experimental simulation of e + e~ — ► X\Xi 
production at y/s = 700 GeV, assuming the integrated luminosity of C = 500 fb" 1 . With beam polar- 
izations (T > e -,V e +) = (—0.8, +0.6) the cross section amounts to a = 150 fb. The event signature is a 
reconstructed hadronic decay W — > qq', a lepton from the decay W — ► Iv [l = e, fx, r) to select clean 
events, and large missing energy of E mms > y/s/2. With a typical selection efficiency e ~ 0.6 and a 
combined branching ratio B = 0.44 prolific event rates are expected. Background from other SM or 
SUSY processes is estimated to be small and will not be considered further. The chargino sample may be 
tripled by including events where both W's are allowed to decay to hadrons. However, the background 
will also increase due to false combinations of jets in reconstructing the two W's and due to background 
from X2X2 production (see above). Such a study goes beyond the aim of the present paper. 

The basis of the analysis is Eq. (4.38) which relates the W energy E\y in the laboratory system with 
the decay angle cos 0^± in the xf res t frame. 

Figures 14(a) and (b) show the angular distribution and energy spectrum for the hypothetical case 
that no QED radiation degrades the nominal production energy. The linear cos 0^± dependence is 
clearly seen at generator level as well as after detector simulation. A fit of the data to the function 
d<j/dcos0* w± ~ H-o?cosfl^±+a§cos2^ ± yields a? = -0.213 ±0.010 and a% = -0.001 ±0.010. This 
value is consistent with the theoretical expectation of Eq. (4.39) with (Kyy±) = —0.216 and demonstrates 
that distortions due to event selection criteria and detector effects are small. The same tendency is 
observed in the energy distribution Fig. 14(b) which falls linearly with Eyy, as compared with a flat 
distribution for unpolarized charginos. 

In the more realistic situation that initial state photon radiation (ISR) and beamstrahlung decrease 
the X\X\ production energy, the angular distribution is no longer linearly falling, as shown in Fig. 14(c). 
Considerable depletions at cos#j^ ± — > ±1 are observed since the constraint E^± is not always valid. 
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Figure 14: Chargino production e + e~ — > XiXi W + X? W~x± with decay angle distributions 
cos in the X\ f^st frame (left panels) and energy distributions Ey/ in the laboratory system (right 
panels). The plots (a) and (b) do not include QED radiation effects while the plots (c) and (d) include 
initial state photon radiation and beamstrahlung. Spectra at generator level (full curves and histograms) 
are compared with events simulated and reconstructed in the detector (dots), the dashed curves indicate 
fits to the data. The simulation is based on polarized beams with (7 ? e -,7 y e +) = (—0.8, +0.6) at ^/s = 
700 GeV and C = 500 fb^ 1 . 



However, since both ISR and beamstrahlung effects can be calculated theoretically and measured pre- 
cisely, they can be unfolded from the data, e.g. by applying a bin-by-bin correction (like in the present 
analysis) or a matrix inversion procedure. Fitting of the QED corrected angular distribution (not shown) 
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to the form 



1 + a\ cos 9w± + a 2 cos 2Q* W± 



(4.40) 



d cos Qt 



results in coefficients 



ai = -0.203 ± 0.020 and a 2 = -O.OOf ± 0.020 



(4.41) 



These values are consistent with the input parameters and confirm with high precision the linear de- 
pendence on cos 6yy characteristic for polarized spin 1/2 chargino production, while higher spin- J states 
would generate the angular distribution 



A sensitivity of a few percents to any term in addition to the linear term can be reached which is an 
important bound in discriminating against higher spin J > 1/2 states. Similar distortions due to QED 
radiation can be seen in the energy spectrum of Fig. 14(d) which is shifted towards lower values and is 
considerably depopulated at the maximum energy. 

(b) Neutralino x%'- 

The distribution of the polar angle in neutralino production e + e~ — ► X2X2 with two identical Majorana 
particles in the final state is forward-backward symmetric. The x.2 polarization, being non-zero for fixed 
polar angle, is asymmetric if the angle is varied from the forward to the backward direction [33]. The 
polarization degree can be enhanced by using polarized electrons/positrons beams. 

The X2 polarization can be determined [45] in the two-body decay X2 ~^ ^Xi if the Z polarization is 
measured in the leptonic decays Z — ► £ + £~. The measurement can also be performed for the hadronic 
decays Z — > cc and bb with c and b flavor tagging. 11 

(Alternatively, if kinematically accessible,) the two-body leptonic decay X2 ~^ i^l^ can provide 
a powerful instrument for determining the x% spin. The x% momenta can be reconstructed, event by 
event, in X2 P & i r production for sequential X2 leptonic decays because the two unknown Xi momenta 
can be fixed by measuring the four visible lepton momenta in the cascade decays x\ ~^ — ► (- + ^~x\ 
and x\ ~^ i'^^ — ► l' + £'~Xi- Furthermore, the slepton mode is a perfect polarization analyzer of 
the decaying neutralino. Explicitly, the angular distribution in the rest frame of the decaying spin-1/2 
neutralino X2 ^ s given by [see Ref. [37]] 



11 A large number of events will be required for the leptonic decays because of the small branching ratio (~ 0.07 for e//x) 
and the small analysis power (~— 0.15) as a result of the almost pure axial-vector Ztt coupling. By contrast, the hadronic 
decays have four times larger branching ratios and much larger analysis powers (~ —0.92 and —0.67 for b and c quarks) 
than the leptonic decays. 
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(4.43) 
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where V^o is the degree of longitudinal x\ polarization and 0| T the angle of the £^ momentum in the 
X2 res t frame with respect to the X2 momentum direction. 

Therefore, the decays x.2 ~~ * ^Xi and/or X2 ~~ *■ ^ d° provide a unique signal for the spin J = 1/2 of 
the neutralino \\. 

Table 2: The threshold behavior and the angular distribution of SUSY and UED particle pair produc- 
tion, and the general characteristics of spin-J field theories. All the characteristics refer to diagonal 
pair production. B and Fq^m generically denote bosons and Dirac, Majorana fermions; [s] s-channel 
exchange only, [s,t, u] potentially all three exchange mechanisms. The parameters k [k 7^ — 1/ depend 
on mass ratios and particle velocities (3; note that, especially, k = 1 for Majorana fermions and for 
self-conjugate bosons [spin > 1/ in [s]-channels - Notice the uniqueness of spin- assignments by mea- 
surements of the polar angle distribution in the slepton sector. Neither threshold excitation nor angular 
distributions are sufficient in the chargino /neutralino sector and final state analyses must be performed 
to determine the spin- 1/2 quantum numbers. 
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5 Summary 

It is apparent from the preceding discussion that the model-independent determination of the spin quan- 
tum numbers of supersymmetric particles is a complex task, with the degree of complexity depending 
on the nature of the particle. Threshold excitation and angular distributions in pair production as well 
as angular correlations in particle decays provide the signals for experimental spin measurements. 
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The predictions for the threshold excitation and the angular distributions in the production processes 
of supersymmetric particles are summarized in Table 2. They are confronted with predictions for 
particles in models of universal extra space dimensions and with general analyses based on the non- 
vanishing of the magnetic dipole moments of all spin > particles. 

Examining these results it turns out that the sin 2 6 law for the production of spin-0 sleptons [for 
selectrons close to threshold] is a unique signal of the spin-0 character. While the observation of the 
sin 2 9 angular distribution is sufficient for sleptons, the /3 3 onset of the excitation curve is a necessary 
but not a sufficient condition for the spin-0 character. Thus the spin determination in the slepton sector 
is conceptually very simple at e + e~ colliders. 

This simple pattern in the slepton sector must be contrasted with the more involved pattern in the 
spin- 1/2 chargino/neutralino sector. Neither the onset of excitation curves nor the angular distributions 
in the production processes provide unique signals of the spin quantum numbers. However, decay angular 
distributions, ~ \d J (9*)\ 2 , do provide a unique signal for the chargino/neutralino spin J = 1/2, albeit 
at the expense of more involved experimental analyses. Using polarized electron/positron beams will in 
general assure that the decaying spin-1/2 particle is polarized; reasonable polarization analysis power 
is guaranteed in many decay processes. 

In toto. The spin of sleptons and charginos/neutralinos can be determined in a model- independent 
way at e + e~ colliders. Similar methods as elaborated for sleptons can be applied in the squark sector 
while gluinos will demand a methodologically separate analysis. 
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Appendix: Cross sections with polarized beams 

Polarized electron and positron beams at e + e~ colliders are useful for diagnosing the properties of 
supersymmetric particles and for unraveling the underlying structure of the SUSY theory [32]. In this 
Appendix we present the general formulae for the production cross sections of R-type smuon/electron 
pairs, and chargino and neutralino pairs in e + e~ annihilation with polarized electron and positron 
beams. 
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For longitudinal beam polarizations V e ± the polarized production cross sections for R-type smuon 
and selectron pairs in e + e~ annihilation are given in terms of the charges Qr and Ql by 

da . ~ l, ~_ 3 



dcosfl WR 



} 



a 3 3 sin 2 9 
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(1 - V e -V e+ )(\Q R \ Z + \Ql\ ) + (TV - 7> e +)(IQ*r - IQlI ) 



(A.l) 



where £ij = £lr, <zr and do = 47ra 2 /3s is the standard normalization cross section of e + e~ annihilation. 
The expressions of the generalized charges Q r and Ql for R-type smuon- and selectron-pair production 
can be found in Eqs. (2.3/2.2) and (3.4/3.3), respectively. For right/left polarized electrons V e - = ± 
and unpolarized positrons, the production cross sections 



dcos9 l R R1 
1 -Mi) 



16 
3 



dcos9 L ' nnJ 16 
project out the bilinear R and L charges Qr and Ql- 



a 3 3 sin 2 0\Q R \ 
do/3 3 sin 2 #|Q L | 



(A.2) 
(A.3) 



The production cross sections for chargino- and neutralino-pairs in e + e annihilation with polarized 
electron and positron beams are given by 



d<T (~ - \ 3 a 
{XiXj} = o°"0P 



dcos 



(1 - V e -V e+ ) {[1 - - fi 2 ) 2 + 3 2 cos 2 9\Q X + 4viVjQ 2 + 2(3cos6Q 3 } 



+ {V e - - V e+ ) {[1 - (m 2 - $Y + IT cos 2 9}Q' l + AfMUjQ^ + 23cos9Q' 3 } 



(A.4) 



where the P-even and P-odd quartic charges Qi and Q\ (i = 1,2,3) are defined in terms of the bilinear 
charges Q a p (a, 3 = L,R) as 

QS = \ [\Qrr\ 2 + \Qrl\ 2 ± \Qlr\ 2 ± \Qll\ 2 } 



-JHe [QrrQrl ± QllQ*lr\ 



Q? } = \ [\Qrr\ 2 - \Qrl\ 2 ± IQl^I 2 t IQllI 2 ] 



(A.5) 



The explicit form of the bilinear charges for the production of the chargino pair x±Xi an d the 
neutralino pairs x?Xj is given in Eqs. (4.3) and (4.4), respectively. Polarized electrons combined with 
unpolarized positrons, 



dcrR 
dcos 9 



dcr L 
dcos 9 



{XiXj} 



{XiXj} 



16 



<j Q 3 



— (Tn 3 

16 



[1 " (^ " V 2 ) 2 + /? 2 cos 2 9] (\Qrr\ 2 + \Q RL \ 2 ) 
+ 8^^ ^(QrrQrl) + 23 cos 9 (\Q rr \ 2 - \Q RL \ 2 ) 
[1 " (tf ~ + P 2 cos 2 9} (\Q L r\ 2 + \Qll\ 2 ) 



+ 8fi tH ^e(Q L RQl L ) + 23 cos 9(\Q L r\ 2 - \Qll\ 2 ) 
project out the bilinear charges Qr^ and Qlu for the x chirality k = R, L. 



(A.6) 



(A.7) 
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